Amino acid sequence homology of monoclonal serum free light chain dimers and tissue deposited light chains in AL amyloidosis: a pilot study.

OBJECTIVES: Diagnosis of light chain amyloidosis (AL) requires demonstration of amyloid deposits in a tissue biopsy followed by appropriate typing. Previous studies demonstrated increased dimerization of monoclonal serum free light chains (FLCs) as a pathological feature of AL. To further examine the pathogenicity of FLC, we aimed at testing amino acid sequence homology between circulating and deposited light chains (LCs). METHODS: Matched tissue biopsy and serum of 10 AL patients were subjected to tissue proteomic amyloid typing and nephelometric FLC assay, respectively. Serum FLC monomers (M) and dimers (D) were analyzed by Western blotting (WB) and mass spectrometry (MS). RESULTS: WB of serum FLCs showed predominance of either κ or λ type, in agreement with the nephelometric assay data. Abnormal FLC M-D patterns typical of AL amyloidosis were demonstrated in 8 AL-λ patients and in one of two AL-κ patients: increased levels of monoclonal FLC dimers, high D/M ratio values of involved FLCs, and high ratios of involved to uninvolved dimeric FLCs. MS of serum FLC dimers showed predominant constant domain sequences, in concordance with the tissue proteomic amyloid typing. Most importantly, variable domain sequence homology between circulating and deposited LC species was demonstrated, mainly in AL-λ cases. CONCLUSIONS: This is the first study to demonstrate homology between circulating FLCs and tissue-deposited LCs in AL-λ amyloidosis. The applied methodology can facilitate studying the pathogenicity of circulating FLC dimers in AL amyloidosis. The study also highlights the potential of FLC monomer and dimer analysis as a non-invasive screening tool for this disease.

Crystal structure of adenylosuccinate lyase from the thermophilic bacterium Thermus thermophilus HB8.

Adenylosuccinate lyase (PurB) catalyzes two distinct reactions in the purine nucleotide biosynthetic pathway using the same active site. The ability to recognize two different sets of substrates is of structural and evolutionary interest. In the present study, the crystal structure of PurB from the thermophilic bacterium Thermus thermophilus HB8 (TtPurB) was determined at a resolution of 2.38 Šby molecular replacement using a structure predicted by AlphaFold2 as a template. The asymmetric unit of the TtPurB crystal contained two TtPurB molecules, and some regions were disordered in the crystal structure. The disordered regions were the substrate-binding site and domain 3. TtPurB forms a homotetramer and the monomer is composed of three domains (domains 1, 2 and 3), which is a typical structure for the aspartase/fumarase superfamily. Molecular dynamics simulations with and without substrate/product were performed using a full-length model of TtPurB which was obtained before deletion of the disordered regions. The substrates and products were bound to the model structures during the MD simulations. The fluctuations of amino-acid residues were greater in the disordered regions and became smaller upon the binding of substrate or product. These results demonstrate that the full-length model obtained using AlphaFold2 can be used to generate the coordinates of disordered regions within the crystal structure.

Identification of a covert evolutionary pathway between two protein folds.

Although homologous protein sequences are expected to adopt similar structures, some amino acid substitutions can interconvert α-helices and ß-sheets. Such fold switching may have occurred over evolutionary history, but supporting evidence has been limited by the: (1) abundance and diversity of sequenced genes, (2) quantity of experimentally determined protein structures, and (3) assumptions underlying the statistical methods used to infer homology. Here, we overcome these barriers by applying multiple statistical methods to a family of ~600,000 bacterial response regulator proteins. We find that their homologous DNA-binding subunits assume divergent structures: helix-turn-helix versus α-helix + ß-sheet (winged helix). Phylogenetic analyses, ancestral sequence reconstruction, and AlphaFold2 models indicate that amino acid substitutions facilitated a switch from helix-turn-helix into winged helix. This structural transformation likely expanded DNA-binding specificity. Our approach uncovers an evolutionary pathway between two protein folds and provides a methodology to identify secondary structure switching in other protein families.

KinFams: De-Novo Classification of Protein Kinases Using CATH Functional Units.

Protein kinases are important targets for treating human disorders, and they are the second most targeted families after G-protein coupled receptors. Several resources provide classification of kinases into evolutionary families (based on sequence homology); however, very few systematically classify functional families (FunFams) comprising evolutionary relatives that share similar functional properties. We have developed the FunFam-MARC (Multidomain ARchitecture-based Clustering) protocol, which uses multi-domain architectures of protein kinases and specificity-determining residues for functional family classification. FunFam-MARC predicts 2210 kinase functional families (KinFams), which have increased functional coherence, in terms of EC annotations, compared to the widely used KinBase classification. Our protocol provides a comprehensive classification for kinase sequences from >10,000 organisms. We associate human KinFams with diseases and drugs and identify 28 druggable human KinFams, i.e., enriched in clinically approved drugs. Since relatives in the same druggable KinFam tend to be structurally conserved, including the drug-binding site, these KinFams may be valuable for shortlisting therapeutic targets. Information on the human KinFams and associated 3D structures from AlphaFold2 are provided via our CATH FTP website and Zenodo. This gives the domain structure representative of each KinFam together with information on any drug compounds available. For 32% of the KinFams, we provide information on highly conserved residue sites that may be associated with specificity.

CATHe: detection of remote homologues for CATH superfamilies using embeddings from protein language models.

MOTIVATION: CATH is a protein domain classification resource that exploits an automated workflow of structure and sequence comparison alongside expert manual curation to construct a hierarchical classification of evolutionary and structural relationships. The aim of this study was to develop algorithms for detecting remote homologues missed by state-of-the-art hidden Markov model (HMM)-based approaches. The method developed (CATHe) combines a neural network with sequence representations obtained from protein language models. It was assessed using a dataset of remote homologues having less than 20% sequence identity to any domain in the training set. RESULTS: The CATHe models trained on 1773 largest and 50 largest CATH superfamilies had an accuracy of 85.6 ± 0.4% and 98.2 ± 0.3%, respectively. As a further test of the power of CATHe to detect more remote homologues missed by HMMs derived from CATH domains, we used a dataset consisting of protein domains that had annotations in Pfam, but not in CATH. By using highly reliable CATHe predictions (expected error rate <0.5%), we were able to provide CATH annotations for 4.62 million Pfam domains. For a subset of these domains from Homo sapiens, we structurally validated 90.86% of the predictions by comparing their corresponding AlphaFold2 structures with structures from the CATH superfamilies to which they were assigned. AVAILABILITY AND IMPLEMENTATION: The code for the developed models is available on https://github.com/vam-sin/CATHe, and the datasets developed in this study can be accessed on https://zenodo.org/record/6327572. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transcriptome Sequencing and Comparison of Venom Glands Revealed Intraspecific Differentiation and Expression Characteristics of Toxin and Defensin Genes in Mesobuthus martensii Populations.

Mesobuthus martensii, a famous and important Traditional Chinese Medicine has a long medical history and unique functions. It is the first scorpion species whose whole genome was sequenced worldwide. In addition, it is the most widespread and infamous poisonous animal in northern China with complex habitats. It possesses several kinds of toxins that can regulate different ion channels and serve as crucial natural drug resources. Extensive and in-depth studies have been performed on the structures and functions of toxins of M. martensii. In this research, we compared the morphology of M. martensii populations from different localities and calculated the COI genetic distance to determine intraspecific variations. Transcriptome sequencing by RNA-sequencing of the venom glands of M. martensii from ten localities and M. eupeus from one locality was analyzed. The results revealed intraspecific variation in the expression of sodium channel toxin genes, potassium channel toxin genes, calcium channel toxin genes, chloride channel toxin genes, and defensin genes that could be related to the habitats in which these populations are distributed, except the genetic relationships. However, it is not the same in different toxin families. M. martensii and M. eupeus exhibit sexual dimorphism under the expression of toxin genes, which also vary in different toxin families. The following order was recorded in the difference of expression of sodium channel toxin genes: interspecific difference; differences among different populations of the same species; differences between sexes in the same population, whereas the order in the difference of expression of potassium channel toxin genes was interspecific difference; differences between both sexes of same populations; differences among the same sex in different populations of the same species. In addition, there existed fewer expressed genes of calcium channel toxins, chloride channel toxins, and defensins (no more than four members in each family), and their expression differences were not distinct. Interestingly, the expression of two calcium channel toxin genes showed a preference for males and certain populations. We found a difference in the expression of sodium channel toxin genes, potassium channel toxin genes, and chloride channel toxin genes between M. martensii and M. eupeus. In most cases, the expression of one member of the toxin gene clusters distributed in series on the genome were close in different populations and genders, and the members of most clusters expressed in same population and gender tended to be the different. Twenty-one toxin genes were found with the MS/MS identification evidence of M. martensii venom. Since scorpions were not subjected to electrical stimulation or other special treatments before conducting the transcriptome extraction experiment, the results suggested the presence of intraspecific variation and sexual dimorphism of toxin components which revealed the expression characteristics of toxin and defensin genes in M. martensii. We believe this study will promote further in-depth research and use of scorpions and their toxin resources, which in turn will be helpful in standardizing the identification and medical applications of Quanxie in traditional Chinese medicine.

DeepRHD: An efficient hybrid feature extraction technique for protein remote homology detection using deep learning strategies.

In computational biology, the Protein Remote homology Detection technique (PRHD) has got undeniable significance. It is mostly important for structure and function identification of a protein sequence. The previous years have seen a challenge that lacks postulating a correlation among the sequences. However, the sequences are of variable length. Thereby, it inhibits the proper derivation of evolutionary information among the sequences. The challenges are the usage of physico-chemical properties as a source to get the evolutionary information and the number of sequences generated every day. This however facilitates a new technique to integrate huge amount of data with a massive feature set. In this article, a new and efficient technique is proposed to predict homology for distantly located sequences of proteins. Deep neural network(CNN-GRU model) is used for the classification of the protein sequences. This is based on different protein families and methods of feature extraction.The efficiency of the proposed model DeepRHD is tested on average 8000 sequences per superfamily taken from SCOP benchmark dataset and the results shows that the proposed model is better than other state of art methods. This model is useful in detecting diseases like sickle cell anemia and influenza and developing a drug thereafter.

Advancing remote homology detection: A step toward understanding and accurately predicting protein function.

Identifying homologous proteins with divergent amino acid sequences can add to our understanding of protein evolution, structure, and function. A new study reports the development of a deep-network-based method to identify 6.8 million new Pfam members, a dramatic singular increase that exceeds a decade of accumulation using traditional approaches.

On the reliability and the limits of inference of amino acid sequence alignments.

MOTIVATION: Alignments are correspondences between sequences. How reliable are alignments of amino acid sequences of proteins, and what inferences about protein relationships can be drawn? Using techniques not previously applied to these questions, by weighting every possible sequence alignment by its posterior probability we derive a formal mathematical expectation, and develop an efficient algorithm for computation of the distance between alternative alignments allowing quantitative comparisons of sequence-based alignments with corresponding reference structure alignments. RESULTS: By analyzing the sequences and structures of 1 million protein domain pairs, we report the variation of the expected distance between sequence-based and structure-based alignments, as a function of (Markov time of) sequence divergence. Our results clearly demarcate the 'daylight', 'twilight' and 'midnight' zones for interpreting residue-residue correspondences from sequence information alone. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome wide identification of MADS box gene family in Musa balbisiana and their divergence during evolution.

MADS box gene family is transcription factor gene family that is involved in growth and development of eukaryotes. In plants the MADS box gene family is mainly associated with floral meristem identity and flower development, apart from being involved in nearly all the phases of plant growth. The MADS box gene family has also been shown to be involved during fruit development and ripening. In this study the MADS box gene family from Musa balbisiana was identified and the divergence of this gene family between Musa balbisiana and Musa acuminata studied. A total of 97 MADS box genes were identified from the genome of Musa balbisiana. Phylogenetic analysis showed that the MbMADS box genes were categorised into type I (α and γ; the ß group was not distinguishable) and type II groups (MIKCc and MIKC* and MIKCc was further divided into 13 subfamilies). The typeII group has the largest number of genes and also showed the most expansion which could be correlated with the whole genome duplications. There were significant differences in the MADS box genes from Musa acuminata and Musa balbisiana during evolution that can be correlated with different floral phenotype and fruit ripening pattern. The divergence of the MADS RIN genes in Musa balbisiana as compared to Musa acuminata might play an important role in the slow ripening of Musa balbisiana fruits.

Cross-Recognition of SARS-CoV-2 B-Cell Epitopes with Other Betacoronavirus Nucleoproteins.

The B and T lymphocytes of the adaptive immune system are important for the control of most viral infections, including COVID-19. Identification of epitopes recognized by these cells is fundamental for understanding how the immune system detects and removes pathogens, and for antiviral vaccine design. Intriguingly, several cross-reactive T lymphocyte epitopes from SARS-CoV-2 with other betacoronaviruses responsible for the common cold have been identified. In addition, antibodies that cross-recognize the spike protein, but not the nucleoprotein (N protein), from different betacoronavirus have also been reported. Using a consensus of eight bioinformatic methods for predicting B-cell epitopes and the collection of experimentally detected epitopes for SARS-CoV and SARS-CoV-2, we identified four surface-exposed, conserved, and hypothetical antigenic regions that are exclusive of the N protein. These regions were analyzed using ELISA assays with two cohorts: SARS-CoV-2 infected patients and pre-COVID-19 samples. Here we describe four epitopes from SARS-CoV-2 N protein that are recognized by the humoral response from multiple individuals infected with COVID-19, and are conserved in other human coronaviruses. Three of these linear surface-exposed sequences and their peptide homologs in SARS-CoV-2 and HCoV-OC43 were also recognized by antibodies from pre-COVID-19 serum samples, indicating cross-reactivity of antibodies against coronavirus N proteins. Different conserved human coronaviruses (HCoVs) cross-reactive B epitopes against SARS-CoV-2 N protein are detected in a significant fraction of individuals not exposed to this pandemic virus. These results have potential clinical implications.

Crystal structures of TTHA1265 and TTHA1264/TTHA1265 complex reveal an intrinsic heterodimeric assembly.

Zinc peptidase M16 family members are widely distributed in most prokaryotic and eukaryotic organisms. M16 family has been divided into three subfamilies, M16A, M16B and M16C, based on sequence alignments and subunit connectivity. TTHA1264, an M16B protein found in Thermus thermophiles HB8, possesses an HXXEH motif essential for Zn2+ binding and catalytic activity. TTHA1265 is another member of M16B, which lacks the metal-binding motif but with a conserved active-site R/Y pair commonly found in the C-terminal half of M16 enzymes. Sequence analysis showed that two genes coding for TTHA1264 and TTHA1265 assemble into a single operon in the bacterial genome. Here, we report the crystal structure of TTHA1265 and TTHA1264/TTHA1265 complex from T. thermophilus HB8. Interestingly, when TTHA1264 and TTHA1265 are present alone, TTHA1264 forms a monomer, TTHA1265 forms a homodimer, respectively. However, TTHA264 and TTHA1265 assembled into a heterodimeric complex, indicating that they prefer to form heterodimer. Biochemical data further confirmed the heterodimeric assembly indicating intrinsic heterodimeric assembly of TTHA1264 and TTHA1265. This property of TTHA1264 and TTHA1265 is consistent with the characteristics of the M16B family.

The ER membrane complex (EMC) can functionally replace the Oxa1 insertase in mitochondria.

Two multisubunit protein complexes for membrane protein insertion were recently identified in the endoplasmic reticulum (ER): the guided entry of tail anchor proteins (GET) complex and ER membrane complex (EMC). The structures of both of their hydrophobic core subunits, which are required for the insertion reaction, revealed an overall similarity to the YidC/Oxa1/Alb3 family members found in bacteria, mitochondria, and chloroplasts. This suggests that these membrane insertion machineries all share a common ancestry. To test whether these ER proteins can functionally replace Oxa1 in yeast mitochondria, we generated strains that express mitochondria-targeted Get2-Get1 and Emc6-Emc3 fusion proteins in Oxa1 deletion mutants. Interestingly, the Emc6-Emc3 fusion was able to complement an Δoxa1 mutant and restored its respiratory competence. The Emc6-Emc3 fusion promoted the insertion of the mitochondrially encoded protein Cox2, as well as of nuclear encoded inner membrane proteins, although was not able to facilitate the assembly of the Atp9 ring. Our observations indicate that protein insertion into the ER is functionally conserved to the insertion mechanism in bacteria and mitochondria and adheres to similar topological principles.

The microprotein Nrs1 rewires the G1/S transcriptional machinery during nitrogen limitation in budding yeast.

Commitment to cell division at the end of G1 phase, termed Start in the budding yeast Saccharomyces cerevisiae, is strongly influenced by nutrient availability. To identify new dominant activators of Start that might operate under different nutrient conditions, we screened a genome-wide ORF overexpression library for genes that bypass a Start arrest caused by absence of the G1 cyclin Cln3 and the transcriptional activator Bck2. We recovered a hypothetical gene YLR053c, renamed NRS1 for Nitrogen-Responsive Start regulator 1, which encodes a poorly characterized 108 amino acid microprotein. Endogenous Nrs1 was nuclear-localized, restricted to poor nitrogen conditions, induced upon TORC1 inhibition, and cell cycle-regulated with a peak at Start. NRS1 interacted genetically with SWI4 and SWI6, which encode subunits of the main G1/S transcription factor complex SBF. Correspondingly, Nrs1 physically interacted with Swi4 and Swi6 and was localized to G1/S promoter DNA. Nrs1 exhibited inherent transactivation activity, and fusion of Nrs1 to the SBF inhibitor Whi5 was sufficient to suppress other Start defects. Nrs1 appears to be a recently evolved microprotein that rewires the G1/S transcriptional machinery under poor nitrogen conditions.

Slipknot or Crystallographic Error: A Computational Analysis of the Plasmodium falciparum DHFR Structural Folds.

The presence of protein structures with atypical folds in the Protein Data Bank (PDB) is rare and may result from naturally occurring knots or crystallographic errors. Proper characterisation of such folds is imperative to understanding the basis of naturally existing knots and correcting crystallographic errors. If left uncorrected, such errors can frustrate downstream experiments that depend on the structures containing them. An atypical fold has been identified in P. falciparum dihydrofolate reductase (PfDHFR) between residues 20-51 (loop 1) and residues 191-205 (loop 2). This enzyme is key to drug discovery efforts in the parasite, necessitating a thorough characterisation of these folds. Using multiple sequence alignments (MSA), a unique insert was identified in loop 1 that exacerbates the appearance of the atypical fold-giving it a slipknot-like topology. However, PfDHFR has not been deposited in the knotted proteins database, and processing its structure failed to identify any knots within its folds. The application of protein homology modelling and molecular dynamics simulations on the DHFR domain of P. falciparum and those of two other organisms (E. coli and M. tuberculosis) that were used as molecular replacement templates in solving the PfDHFR structure revealed plausible unentangled or open conformations of these loops. These results will serve as guides for crystallographic experiments to provide further insights into the atypical folds identified.

Multidimensional Phylogenetic Metrics Identify Class I Aminoacyl-tRNA Synthetase Evolutionary Mosaicity and Inter-Modular Coupling.

The role of aminoacyl-tRNA synthetases (aaRS) in the emergence and evolution of genetic coding poses challenging questions concerning their provenance. We seek evidence about their ancestry from curated structure-based multiple sequence alignments of a structurally invariant "scaffold" shared by all 10 canonical Class I aaRS. Three uncorrelated phylogenetic metrics-mutation frequency, its uniformity, and row-by-row cladistic congruence-imply that the Class I scaffold is a mosaic assembled from successive genetic sources. Metrics for different modules vary in accordance with their presumed functionality. Sequences derived from the ATP- and amino acid- binding sites exhibit specific two-way coupling to those derived from Connecting Peptide 1, a third module whose metrics suggest later acquisition. The data help validate: (i) experimental fragmentations of the canonical Class I structure into three partitions that retain catalytic activities in proportion to their length; and (ii) evidence that the ancestral Class I aaRS gene also encoded a Class II ancestor in frame on the opposite strand. A 46-residue Class I "protozyme" roots the Class I tree prior to the adaptive radiation of the Rossmann dinucleotide binding fold that refined substrate discrimination. Such rooting implies near simultaneous emergence of genetic coding and the origin of the proteome, resolving a conundrum posed by previous inferences that Class I aaRS evolved after the genetic code had been implemented in an RNA world. Further, pinpointing discontinuous enhancements of aaRS fidelity establishes a timeline for the growth of coding from a binary amino acid alphabet.

Synthetic virions reveal fatty acid-coupled adaptive immunogenicity of SARS-CoV-2 spike glycoprotein.

SARS-CoV-2 infection is a major global public health concern with incompletely understood pathogenesis. The SARS-CoV-2 spike (S) glycoprotein comprises a highly conserved free fatty acid binding pocket (FABP) with unknown function and evolutionary selection advantage1,2. Deciphering FABP impact on COVID-19 progression is challenged by the heterogenous nature and large molecular variability of live virus. Here we create synthetic minimal virions (MiniVs) of wild-type and mutant SARS-CoV-2 with precise molecular composition and programmable complexity by bottom-up assembly. MiniV-based systematic assessment of S free fatty acid (FFA) binding reveals that FABP functions as an allosteric regulatory site enabling adaptation of SARS-CoV-2 immunogenicity to inflammation states via binding of pro-inflammatory FFAs. This is achieved by regulation of the S open-to-close equilibrium and the exposure of both, the receptor binding domain (RBD) and the SARS-CoV-2 RGD motif that is responsible for integrin co-receptor engagement. We find that the FDA-approved drugs vitamin K and dexamethasone modulate S-based cell binding in an FABP-like manner. In inflammatory FFA environments, neutralizing immunoglobulins from human convalescent COVID-19 donors lose neutralization activity. Empowered by our MiniV technology, we suggest a conserved mechanism by which SARS-CoV-2 dynamically couples its immunogenicity to the host immune response.

Mass Screening of SARS-CoV-2 Variants using Sanger Sequencing Strategy in Hiroshima, Japan.

This study aimed to develop the feasible and effective universal screening strategy of the notable SARS-CoV-2 variants by Sanger Sequencing Strategy and then practically applied it for mass screening in Hiroshima, Japan. A total of 734 samples from COVID-19 confirmed cases in Hiroshima were screened for the notable SARS-CoV-2 variants (B.1.1.7, B.1.351, P.1, B.1.617.2, B.1.617.1, C.37, B.1.1.529, etc.). The targeted spike region is amplified by nested RT-PCR using in-house designed primer set hCoV-Spike-A and standard amplification protocol. Additionally, randomly selected 96 samples were also amplified using primer sets hCoV-Spike-B and hCoV-Spike-C. The negative amplified samples were repeated for second attempt of amplification by volume-up protocol. Thereafter, the amplified products were assigned for Sanger sequencing using corresponding primers. The positive amplification rate of primer set hCoV-Spike-A, hCoV-Spike-B and hCoV-Spike-C were 87.3%, 83.3% and 93.8% respectively for standard protocol and increased to 99.6%, 95.8% and 96.9% after second attempt by volume-up protocol. The readiness of genome sequences was 96.9%, 100% and 100% respectively. Among 48 mutant isolates, 26 were B.1.1.7 (Alpha), 7 were E484K single mutation and the rest were other types of mutation. Moreover, 5 cluster cases with single mutation at N501S were firstly reported in Hiroshima. This study indicates the reliability and effectiveness of Sanger sequencing to screen large number of samples for the notable SARS-CoV-2 variants. Compared to the Next Generation Sequencing (NGS), our method introduces the feasible, universally applicable, and practically useful tool for identification of the emerging variants with less expensive and time consuming especially in those countries where the NGS is not practically available. Our method allows not only to identify the pre-existing variants but also to examine other rare type of mutation or newly emerged variants and is crucial for prevention and control of pandemic.

SIRT2 regulates proliferation and chemotherapy response of MLL-ENL-driven acute myeloid leukemia.

Both MLL-AF9 and MLL-ENL leukemia fusion proteins drive oncogenic transformation of hematopoietic cells through their N-terminal DNA/histone binding mixed-lineage leukemia 1 domain and C-terminal fragment of AF9 or ENL containing an unstructured linker region and the ANC1 homology domain, which recruits transcription factors. Despite of their structural similarity, acute myeloid leukemia (AML) patients bearing MLL-ENL show more adverse outcomes compared to those with MLL-AF9. We recapitulated the clinical patterns of these two MLL-fusions driven AMLs using murine models and found that MLL-ENL AML cells showed slower cell cycle progression and more resistance to standard chemotherapy than MLL-AF9 cells. These phenotypes were primarily controlled by the linker regions of ENL and a highly conserved lysine residue K469 within. Substitution of K469 with an acetylated mimic glutamine abolished the ability of MLL-ENL to suppress proliferation and promote chemo-resistance. We showed that deacetylase Sirt2 might act as an upstream regulator of MLL-ENL. Deletion of Sirt2 promoted proliferation of AML cells with either MLL fusions. Importantly, loss of Sirt2 greatly enhanced the sensitivity of the MLL-ENL AML cells to chemo-treatment. Taken together, our study uncovered a unique regulatory role of Sirt2 in leukemogenesis and suggested targeting SIRT2 as a new way to sensitize MLL-ENL AML patience for chemotherapy.

Application of Homology Modeling by Enhanced Profile-Profile Alignment and Flexible-Fitting Simulation to Cryo-EM Based Structure Determination.

Application of cryo-electron microscopy (cryo-EM) is crucially important for ascertaining the atomic structure of large biomolecules such as ribosomes and protein complexes in membranes. Advances in cryo-EM technology and software have made it possible to obtain data with near-atomic resolution, but the method is still often capable of producing only a density map with up to medium resolution, either partially or entirely. Therefore, bridging the gap separating the density map and the atomic model is necessary. Herein, we propose a methodology for constructing atomic structure models based on cryo-EM maps with low-to-medium resolution. The method is a combination of sensitive and accurate homology modeling using our profile-profile alignment method with a flexible-fitting method using molecular dynamics simulation. As described herein, this study used benchmark applications to evaluate the model constructions of human two-pore channel 2 (one target protein in CASP13 with its structure determined using cryo-EM data) and the overall structure of Enterococcus hirae V-ATPase complex.